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The phenomenon of frustration, which gives 
rise to many fascinating phenomena, is conven- 
tionally associated with the topology of non- 
bipartite lattices, where nearest-neighbor (nn) in- 
teractions and global connectivity compete in the 
lowering of energy. The issue of rattling atoms 
in spacious lattice sites is a separate occurrence 
that can also lead to a high density of low en- 
ergy states (unusual low temperature thermody- 
namics) and to practical applications such as in 
improved thermoelectric materials. In this letter 
we address a unique situation where both phe- 
nomena arise: a four-fold single-site instability 
leads to rattling of cations on a diamond structure 
sublattice where nn interactions frustrate simple 
ordering of the displacements. The system deals 
with this coupling of rattling-|-frustration by com- 
mensurate ordering. Such a disorder-order tran- 
sition may account for the second phase tran- 
sition seen in KOS2O6 within the superconduct- 
ing state, and the unusual low-energy dynamics 
and associated electron-phonon coupling can ac- 
count for the qualitative differences in physical 
properties of KOS2O6 compared to RbOs206 and 
CsOs206, all of which have essentially identical 
average crystal and electronic structures. 

The pyrochlore-lattice-based structure with a potential 
to support magnetic frustration has attracted attention 
to AOS2O6 (A=K, Rb,Cs) group. Unexpectedly large 
variation of the superconducting Tc throughout the group 
(from 3.3 K in CsOs206 to 9.7 K in KOS2O6) 0,1113 to- 
other with reports of anomalous nuclear spin relaxation 
and indications of anisotropic order parameter ^J in 
KOS2O6 pointed to a possibility of unconventional pair- 
ing and fueled the early experimental interest. While the 
issue of superconductivity remains controversial in the 
light of recent pressure experiments Ij, unusual trans- 
port and thermodynamic properties were found in the 
normal state of KOS2O6 in sharp contrast to the standard 
metallic behavior of RbOs206 and CsOs206 1^,13,1^' 

Uniquely to KOS2O6 within this class, the normal-state 
conductivity exhibits a non-Fermi-liquid behavior char- 
acterized by a concave temperature dependency down to 
low temperatures jsl US^- The low temperature linear 
specific heat coefficient is estimated to be substantially 



larger than in RbOs206 and CSOS2O6 |ll|. Recently 
an intriguing A-shaped peak in the specific heat was ob- 
served in good quality KOS2O6 single-crystals indicative 
of a phase transition at Tp = 7 K [lOj, within the su- 
perconducting state. This observation was recently con- 
firmed _9J. Notably, the peak position and shape do not 
change even when the superconductivity is suppressed 
below 7 K by the external field. Insensitivity to such 
a profound change of the electronic state indicates that 
the peak is rooted in the lattice dynamics rather than 
intrinsic electronic degrees of freedom. 

Electronic structure investigations [ij |l2| have re- 
vealed a considerable bandwidth of the Os-5d-t2g 12-band 
complex of about 3 eV which does not support the idea 
of local moment formation on the Os sites nor any emer- 
gence of frustration due to the pyrochlore topology of the 
Os sublattice, made of a three-dimensional network of 
vertex-sharing tetrahedra. Instead we find that a signifi- 
cant frustration, not magnetic but structural, takes place 
on the diamond sublattice occupied by K ions. We have 
shown previously |ll| that the symmetric (A^) potas- 
sium phonon mode is unstable and that the energy can 
be lowered by several meV/atom(K) through rather large 
displacements of the K ions. Here we construct, based 
on first principles calculations, the effective potential de- 
scribing fourfold symmetric displacements of K ions off 
their ideal diamond-lattice sites, with nn coupling leading 
to a highly frustrated system of displacements. Dynami- 
cal simulations for finite clusters reveal a classical ground 
state with complex pattern of displacements. 

In Fig. ^we show the AOS2O6 lattice which consists of 
Os on a pyrochlore sublattice, having one O atom bridg- 
ing each Os nn pair. The cavities in the Os-0 network 
are filled with alkali ions, which themselves form a dia- 
mond lattice, composed of two fee sublattices. Using 
a full-potential linearized augmented-plane-waves code 
Wien2k 13] we have performed a series of calculations 
in which K ions move along the (111) direction: (i) the 
two fee sublattices are displaced in opposite direction 
(symmetric Ag mode), (ii) same as (i) with the O po- 
sitions allowed to relax, (iii) only one fee sublattice is 
displaced. Comparing results (i) and (ii) reveals a non- 
negligible O relaxation only for large K displacement ^j 
(the energy vs displacement curve has slightly less steep 
walls when O ions are allowed to relax). Since the O re- 




weak coupling of about -2 K for each pair of orbitals ac- 
counts for the singlet-triplet splitting. 

Ab initio calculations [ll| revealed much less anhar- 
monicity in the Rb and especially Cs potentials, which 
can treated as a perturbation and neglected for the 
present purposes. Within this approximation Rb (Cs) 
dynamics is described by 3D oscillator with a frequency 
of 44 K (61 K). Recent analysis of specific heat data 
by Briihwiler et al. ;^j led to a frequency of 60 K for 
RbOs206, which we find a reasonable agreement given 
the approximations. Localized modes were previously 
observed in specific heat of RbOs2 06 and CsOs206 by 
Hiroi et al. [l5( The harmonic oscillator root mean square 
displacement is 



FIG. 1: The AOS2O6 lattice. The atomic species are maker 
with different colors: Os (gray), O (red) and alkali metal 
(blue). The pyrochlore sublattice of Os atoms is highlighted. 
Notice the alkali sublattice with diamond structure. 



laxation effect is minor we consider the Os-0 network to 
be rigid for the following discussion. The inter-ionic dis- 
tances together with the geometry of the Os-0 network 
with spacious channels along the nn K-K bonds suggest 
nn coupling to dominate over longer range interaction. 
The effective Hamiltonian becomes 
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where the first term is the on-site Hamiltonian and sec- 
ond describes the nn coupling (interaction). The on-site 
potential, which captures the essential tetrahedral local 
symmetry, consists of a spherical and the next non-zero 
term in spherical harmonic expansion (J^jm)) while the ra- 
dial dependency is described by even and odd 6th-order 
polynomials Pei^i) and PoUi) obtained by fitting the ab 
initio data from type (iii) calculations (Fig. |2J). 

We have solved the quantum-mechanical single site 
problem by numerical integration on real space grid (de- 
tails of the calculation can be found in Ref. Il4j) . The 
low energy spectrum up to 80 K (containing 20 states) is 
characterized by a singlet-triplet split ground state (8 K 
splitting) separated by a gap of about 25 K from excited 
states. This essential difference from the harmonic poten- 
tial with singlet ground state is reflected also by a Schot- 
tkyesque anomaly in the single-site specific heat.[lj| The 
sharpness of the observed peak, however, rules out the 
Schottky anomaly scenario, pointing instead to a collec- 
tive transformation involving inter-site coupling. A par- 
ticularly useful way of looking at the quasi-degenerate 
quadruplet ground state, is in terms of four symmetry 
related local orbitals centered in the local minima. A 
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and we obtain Afu, = 0.1 5ao (Bohr radius) and Acs — 
O.loo for the mean square displacement at zero temper- 
ature, which we will use below in estimation of strength 
of the inter-site coupling. 

The interaction can be obtained by following the force 
acting on a fixed ion when its neighbors are uniformly 
displaced. The simplest form of central pair force that 
describes reasonably well the ab initio data is F{r) = 
A- + Br, corresponding to a pair potential T^(ri,r2) — 



A\ri2 



|ri2p (rj are ion coordinates) . In Fig. [2Iinset) 



we show the first principles force together with the model 
fit, the values A=-88 mRy/ao and B=7.9 mRy/ap are es- 
sentially the same for all three oxides. As expected from 
its electrostatic origin the pair force is repulsive for ad- 
missible values of r. Using the electrostatic force — ;^ 
instead of an ad hoc Taylor expansion, the force in Fig. 
m would not crossover to the positive values, but only 
reach zero for zero displacement. The observed behav- 
ior is qualitatively consistent with faster decay of the 
interaction due to screening. The interaction Wij{^i,^j) 
between ions at sites Ri and Rj is obtained after sub- 
traction of contributions accounted for in the on-site po- 
tential: 



W^j (li, ^j) -l^(Ri + ii, Rj + ^j)- ^(Ri, Rj 



^j) 
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The result is a directional (non-central) potential, whose 
dipolar form becomes clear in the small displacement 
limit: 



m,i^i,^j)^A 



(Rjj • $i)(Rij ■ $.j] 



-(A_ 



-B)$i-$j. (4) 



We start the discussion of nn coupling in the sim- 
pler quasiharmonic case (RbOs206 and CsOs206). Small 
mean displacements justify the use of a dipolar approxi- 
mation Q to estimate the interaction energy. Using the 
Schwarz inequality |(^i • ^j)P < {£,i){£,'j) the interaction 
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FIG. 2: The on-site potential. The on-site potential along 
the nn bond direction (the nn site is in the direction of positive 
X-axis). The vertical line denote the eigenenergies on the same 
scale, the thickness represents degeneracy ranging from 1 to 
3. In the inset force acting on undisplaced ion as a function 
of the uniform displacement of its neighbors is shown. The 
symbol represent the ab initio data, the full line is the fit with 
linear pair force. 



energy per site for a two-site problem has the following 
upper bound 
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yielding 23 K and 10 K for Rb and Cs respectively. Going 
one step further and considering the problem of four sites 
connected to their common neighbor the upper bound is 
reduced for purely geometrical reasons to 16 K and 7 
K respectively. Smallness of the interaction energy in 
comparison to the Einstein frequencies and the singlet 
character of the on-site groundstate lead to conclusion 
that Rb and Cs dynamics is essentially local. 

KOS2O6 presents the opposite limit of extreme anhar- 
monicity, leading to quasi-degeneracy of the ground state 
and large spatial fluctuations, together with important 
inter-site coupling. In the following we evaluate the ma- 
trix elements of W in the basis of products \a) of the 
on-site eigenstates. The basis of local orbitals for the 
ground state quadruplet has an advantage of keeping the 
Hamiltonian in quasi-diagonal form by maximizing the 
diagonal terms (a/3|W^|/3Q;) and minimizing the leading 
off-diagonal contributions {af3\W\ja) . Moreover there is 
a natural one-to-one mapping between the local orbitals 
and nn bonds, namely we denote an orbital and a bond 
with the same index if the orbital represents displace- 
ment in the direction of the bond. Due to the symme- 
try there are only four independent diagonal matrix ele- 
ments, which can be expressed in the form (in kelvins): 
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where both the bond index R and the orbital indices a, 
(3 run from 1 to 4. These numbers can be understood in 
terms of the approximate formula Q, taking into account 
the inversion symmetry about the bond center. The rel- 
evant off-diagonal terms yield about 15 K (additional 2 
K comes from the on-site Hamiltonian) . The off-diagonal 
terms also provide coupling to products including excited 
states with the largest ones being about 1/3 of the cor- 
responding energy difference, providing thus small but 
non-negligible quantum mechanical coupling. 

Origin of frustration. Building a lattice model from 
bonds © helps to understand frustrated nature of the 
present system. Although it is not justifiable to neglect 
the excited states completely, they will only renormalize 
the parameters without changing the form of the four- 
state form of the Hamiltonian in the low energy sector. 
Building a lattice Hamiltonian from the bonds 10, us- 
ing the fact that the third term yields a constant when 
summed over the bonds, we get an expression 
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The first term of Q is the classical Potts Hamiltonian 
|l6|. (5^a yields 1 when the neighboring sites i, j are oc- 
cupied by the same state and zero otherwise. In the 
second term R{ij) is an index of the bond between sites 
i and j, (5^rcj ■■. yields 1 if orbital on site i corresponds to 
displacement in the direction of the bond ij and zero oth- 
erwise. The bare values of parameters a and b are -162 K 
and 371 K respectively. Since the second term describes 
states whose energy is above the already neglected ex- 
cited states it is consistent to rule these states out by 
putting b equal to -l-oo. The second term thus becomes 
a constraint on admissible configurations and introduces 
frustration into the system. The leading quantum me- 
chanical correction H on-site, bare value of which is an 
order of magnitude smaller than a, is provided by tun- 
neling between the local orbitals. Filling the lattice such 
that we minimize the contribution of an arbitrary first 
site (only 3 bonds can yield a due to the constraint) one 
can readily see that an arrangement with the same en- 
ergy cannot be placed on the neighboring sites. Unlike in 
the case of geometrical frustration of nn antiferromagnets 
no odd-length loops are necessary to produce frustration. 
In fact the above mechanism would apply even to Bethe 
lattice with no loops at all. While we cannot make con- 
clusions about the degeneracy of the groundstate, the 
frustrating constraint is expected to reduce the transition 
temperature below the energy scale defined by parameter 
a. 

Dynamical simulations. While the effective Hamilto- 
nian |(7J| can be useful for investigating general features 



of the phase transition and is well suited for analytical 
approach, in the rest of this paper we pursue a separate, 
purely numerical approach to probe aspects of the order- 
ing that we anticipate at Tp. Addressing this question in 
full generality is very difficult. Insight can be gained by 
minimizing the potential energy, i.e. pursuing the classi- 
cal (large mass M) limit, for finite clusters with periodic 
boundary conditions. This is still a formidable compu- 
tational task due to a large number of local minima. To 
approach and possibly reach the global minimum we have 
used a damped molecular dynamics combined with sim- 
ulated annealing. In particular we have integrated the 
classical equation of motion 
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where J^ is the actual force, the /3(T) c>c Vt is a friction 
parameter and Q(T) is a Gaussian random vector with 
half-width proportional to T. The effective temperature 
T was successively reduced, Ti = eTi_i(e < 1), until 
minimum was reached. 

The minimum of the 1x1x1 (single primitive cell) clus- 
ter can be described as parallel displacement of all ions 
along one of the bond directions with different displace- 
ment values (1.54ao toward the nn site and 1.24ao away 
from nn site) on the two sublattices (the global minimum 
is of course degenerate with respect to the sublattice ex- 
change) . The ordering on a 2 x 2 x 2 cluster is character- 
ized by uniform displacements along different bonds as 
shown in Fig. El The minima for 3x3x3 and larger 
clusters are difficult to understand in real space since the 
displacements are neither uniform nor limited to bond di- 
rections. Nevertheless, common features include a small 
net displacement per sublattice (less then O.lag) and an 
average displacement of 2.0ao per site with a standard de- 
viation of about 0.3ao. The size of displacements is likely 
to be overestimated due to neglect of the kinetic energy, 
effect of which can be qualitatively visualized as replac- 
ing the point particles with probability density clouds. 
Moreover Fourier transform of the displacement vectors 

^Q(R-i) 




Sa{q) = J^Yl exp(«q • Ri)^a(Rz), 



(9) 



FIG. 3: Minimum energy configuration for 2x2x2 cluster. 
Different colors correspond the two fee subfattices, green lines 
mark the nn bonds. (Note that the same dispfacements for 
sites on the same edge of the cube are not enforced by the 
boundary conditions.) 



becomes very flat while the ordering wavevectors are sen- 
sitive to boundary conditions. The absolute value of this 
energy has no relevance for low energy scale of the order- 
ing transition, but its convergence indicates that mini- 
mum energy is being reached. 

Our results provide a picture of potassium dynamics 
in KOS2O6 governed by an effective Hamiltonian charac- 
terized by an unusually soft and broad local potential, 
which allows for large excursions of K ions resulting in 
significant and frustrating nn coupling. To address the 
question of ordering tendencies we have used classical 
simulations for finite clusters which provide a compli- 
cated but distinct pattern with multiple-q ordering and 
large displacements. Since the purely numerical model 
is not well suited for addressing general questions con- 
cerning the phase transition and for understanding the 
essence of the present physics we have also proposed an 
analytic model with only a few parameters. This model 
is formally a three-dimensional ferromagnetic four-state 
Potts model with an additional constraint on possible 
configurations. While the unconstrained model is known 



revealed that there are only a few non-vanishing q- 
components for each cluster size. Even with the lowest 
cooling rate we were not able to obtain the minimum for 
5x5x5 cluster unambiguously, which strongly suggests 
that periodicity of 5 unit cells is not commensurate with 
the ordering tendencies in the system. The results are 
summarized in TableQ] Fourier transforms are character- 
ized by 2 or 3 dominant components with (2/3,2/3,0) and 
(1/2,1/2,1/2) appearing whenever allowed by the cluster 
size. Comparison of the minimum energies for different 
clusters indicates that beyond 3x3x3 the energetics 



TABLE I: Minimum potential energy and the dominant 
Fourier components (only one member of ±q pair is shown) 
for the ground states of clusters of different size (in brackets) . 
The vectors are in the units of — where a=10.101 A. 



cluster size 



1 X 1 X 1 (2) 
2x2x2 (16) 
3x3x3 (54) 
4x4x4 (128) 
6x6x6 (432) 



E (mRy) 



-1.99 
-17.11 
-18.01 
-18.11 
-18.20 
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to exhibit a first order mean-field-like phase transition 
|l6| , the constraint cannot be relieved in a simple way by 
the system and is likely to change behavior of the model. 

Our calculations suggest a natural explanation for the 
second peak observed in the specific heat of KOS2O6 [ifl 
as a phase transition of the potassium sublattice to super- 
cell order. Anomalies of low temperature electronic prop- 
erties such as non- Fermi-liquid conductivity and large lin- 
ear specific heat coefficient |l0j can be explained as con- 
sequence of atomic motion, which does not freeze down 
to the ordering transition at 7 K. We point out that large 
excursions of K ion should affect the NMR measurements 
due to quadrupolar interaction and might be responsible 
for observed anomalies [^] . The dynamics of Rb and Cs 
ions is very different: the local behavior is quite different 
due to the larger ionic radii which gives a different energy 
scale, and this distinction in turn negates the inter-site 
interaction, leaving a simple quasiharmonic local mode. 
This result of our first principle calculations fits well with 
the observed specific heat and conductivity jfll. I1!tI|. 

If it proves possible, synthesis of Ka;Rbi_a;Os206 will 
provide a means of introducing 'vacancies' into the model 
Hamiltonian |(7J. The combination of rattling and frus- 
trating nn interaction, facilitated by 'fine tuning' of the 
potassium ionic radius to the size of osmium-oxygen cage, 
provides a novel physical system, which exhibits a phase 
transition as low temperature. 
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